Prethermalization and dynamical transition in an isolated trapped ion spin chain 
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We propose an experimental scheme to observe prethermalization and dynamical transition in one- 
dimensional XY spin chain with long range interaction and inhomogeneous lattice spacing, which 
can be readily implemented with the recently developed trapped-ion quantum simulator. Local 
physical observables are found to relax to prethermal values at intermediate time scale, followed by 
complete relaxation to thermal values at much longer time. The physical origin of prethermalization 
is explained by spotting a non-trivial structure in lower half of the energy spectrum. The dynamical 
behavior of the system is shown to cross difference phases when the interaction range is continuously 
tuned, indicating the existence of dynamical phase transition. 



The dynamical properties of isolated quantum many- 
body systems have been under intense interest in recent 
years [TJ [5]. On the theory side, the research has been 
centered on whether and how an isolated quantum sys- 
tem approaches thermal equilibrium. While certain ob- 
servables are found to relax to equilibrium in some large 
systems |3HZ], it remains unclear on what conditions and 
time scale equilibration occurs in generic systems |S]- 
ITT] . On the experimental side, recent progress with cold 
atoms |12IfT4"] and trapped ions |15fBU] has made it pos- 
sible to simulate well controlled simple models, such as 
one-dimensional (ID) Bose gas and transverse field Ising 
model. These quantum systems can be well isolated from 
the environmental bath and have long coherence time, 
while their physical properties can be measured at indi- 
vidual atomic level, providing an unprecedented oppor- 
tunity for studying non-equilibrium dynamics in closed 
interacting systems. 

A particularly intriguing phenomenon in this context 
is called prethermalization |21| . which has been shown to 
emerge in various theoretical setups [22-24], and experi- 
mentally observed in cold atomic gas [14]. The emergence 
of prethermalization is characterized by establishment of 
quasi-stationary state at intermediate time scale, and fol- 
lowed by relaxation to stationary state at much longer 
time scale (thermalization) . Physical origin of prether- 
malization, however, is still elusive, and is primarily spec- 
ulated to be related to quasi-integrability of the model 

[HE3]. 

In this paper, we propose a new experimental scheme 
for observing and studying prethermalization and related 
dynamical transition in a XY spin model, which can be 
implemented with the current trapped-ion quantum sim- 
ulator [19 . Our model features long range spin-spin in- 
teraction with inhomogeneous lattice spacing, and unlike 
many other systems, the prethermalization can occur al- 
ready for as few as a dozen of spins, allowing for its ob- 
servation in current experimental systems. The prether- 
malization shown up in this system has a quite different 
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Figure 1: Schematic of the proposed experimental setup: A 
chain of N ions are trapped along the z direction in a ID 
harmonic linear Paul trap. The global Raman beams generate 
spin-dependent force along x direction, resulting in effective 
Ising-type interaction. To induce dynamics, a focused laser 
beam is applied on one end of the ion chain to selectively flip 
only the first spin. 

mechanism and we find that a non-trivial structure in 
the energy spectrum resulting from long-range interac- 
tion and inhomogeneous lattice is responsible for the oc- 
currence of prethermalization in our model. In addition, 
by tuning the range of interaction with an experimental 
knob, we find the dynamical behavior of system exhibits 
three different phases: thermalization only, prethermal- 
ization followed by thermalization, and prethermaliza- 
tion only. The transition between different phases be- 
comes sharper and sharper with increased number of 
spins, hinting the existence of dynamical phase transi- 
tions |25j . 

Model and its dynamics Our spin model is based on 
the experimental system of a chain of ions confined in a 
linear Paul trap (Fig. [I]). Through proper configuration 
of Raman beams, the optical dipole force can generate 
an effective transverse field Ising model [TBI H~7l 125] : 

N N 

H = J2Ji,jO-fv*+Bj2<r?, (1) 

i<j i=l 

where <7j is the spin- 1/2 Pauli matrix for the i th ion qubit. 
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The interaction coefficient Jjj in Eq. (1) is given by 
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where /i is the Raman beatnote frequency, f2 is the effec- 
tive Rabi frequency which is assumed to be uniform on 
all the ions, {w m } are the phonon mode frequencies of 
ions in x direction, and r) i m are the Lamb-Dicke param- 
eters measuring the coupling between the «th ion and the 
mth phonon mode. We are interested in the region where 
B 3> max{ Jij}. In this limit, the o-fo~~j and c~cr~ terms 
in Eq. [TJwill be energetically forbidden, and we end up 
with the XY Hamiltonian: 

H « H XY = ]T 2J itj (a+aj + h.c.) + B^of (2) 

i<j i 

A unique feature of the Hamiltonian ([T| & [2]) realized 
with the ion system is that the interaction characterized 
by Jij is long-ranged and the range of interaction can be 
readily tuned by changing the beatnote frequency /i. In 
particular, in a range of /i, Jij can be roughly approxi- 
mated by an power-law decay with Jij ~ \i — j\~ a , where 
a varies from to 3 when we tune /x[19J. In our follow- 
ing analysis, for a given fx, we fit the coefficient Jij with 
J.^ ~ \i — j\~ a and use the fitting parameter a as an 
indicator of the range of iteration. 

To study dynamics of the model Hamiltonian (2), we 
first initialize all the spins through optical pumping to 
the spin down state with erf = —1, which is an eigenstate 
of Hxy and hence stationary. We then use a focused 
laser beam to flip the first spin (left end ion) to erf = 1 
(see Fig. [TJ. The starting state \ip(0)) = | f44 ' ' ' i) is 
no longer an eigenstate of Hxy and subject to dynamics 
with \tp(t)) = e iHxYt / h \i/)(Q)}. We consider time evolu- 
tion of local observables (erf (t)} and their correlations 
(erf (t)(Tj (t)} which can be directly measured in experi- 
ments. For convenience of description of the dynamics, 
we introduce the operator 
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where the coefficient /j = (i — ^i^-)/(N — 1) is equally 
distributed between [—1,1] from i = 1 to i = N. The ex- 
pectation value of C varies between [—1,1] and physically 
measures the relative position of the spin excitation. It's 
easy to check that (V'(O)|C|'0(O)) = —1, meaning the spin 
excitation is at the left edge of the chain. For any state 
with spatial inversion symmetry around the center of the 
chain, (C) = 0. 

Prethermalization and dynamical transition: To find 
out the dynamical behavior, we first perform numerical 
calculation with an TV — 16 ion chain, which corresponds 
to the size of the current experimental platform for ion 
quantum simulator |19| . As shown in Fig. [2] we pick 
up two parameter settings with the corresponding fitting 




Figure 2: (a) Distribution of Jij for short range interaction 
with the beatnote frequency set at [i — 5.2 MHz, the trap 
frequency in z direction uu z = 100 KHz, and ri x Q — 40 KHz. 
The corresponding fitting parameter a w 2.6 in this case, (b) 
Distribution of Jij for long range interaction with fj, — 5.02 
MHz, uiz = 600KHz, r) x Q. = 3.9KHz, and the corresponding 
fitting parameter a w 0.52. In both cases we have trap fre- 
quency in x direction uj x — 5 MHz and average interaction 
strength J = J t] /N 2 = 20 Hz. 
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Figure 3: (a-b) Short time dynamics of of and C for short- 
range (a) and long-range (b) interaction, (c-d) Long time 
dynamics of time-averaged C and ofcrj for (c) short-range 
and (d) long-range interaction. The parameters for short- 
range and long-range interaction are given by Fig. 2 



parameter a ~ 2.6 and a w 0.52, which represent respec- 
tively short-range and long-range interaction. The distri- 
butions of the exact coupling coefficients J,^ are shown 
in Fig. [2] for these two cases. 

For these choices of parameters, the short time dynam- 
ics with t £ [0,2/ Jq] for all (erf) and (C) are shown in 
Fig. [3] In the short-range interaction case, one observes 
that the spin excitation, initially located at the left edge 
of the chain, almost coherently travels to the other side 
and oscillates back and forth with relatively small dis- 
persion. In contrast, the spin excitation diffuses to the 
rest of the chain in the long-range interaction case and 
somehow get locked before it reaches the middle of the 
chain (with (C) « -0.4). 

To better show the long time dynamics, we use the 
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Figure 4: A dynamical "phase diagram" with regard to the 
interaction range characterized by the fitting parameter a. 
The TV = 16 case uses the parameters specified in Fig. [l]and 
N — 64, 256 cases use the same parameters except that lu z is 
scaled down by lo z oc VmJV /N to maintain chain stability. 



finite-time-averaged quantity A(t). defined as A(t) = 
| J*{A(T))dT for the observable A [571125]. This will av- 
erage out temporal fluctuations on short time scale (the 
following dynamical behaviors are qualitatively the same 
even without performing any time averaging). The long 
time dynamics is shown in Fig. [3]jc-d). In the short- 
range interaction case, the spin excitation position C, as 
well as the spin correlation crfcr?, relax to the stationary 
values at around T rs 10/ Jo = 500ms. In the long-range 
interaction case, observables first reach quasi-stationary 
(prethermal) values at time scale To, and further relax to 
the stationary (thermal) values at much longer time scale 
(10 4 T ). The emergence of prethermalization is mani- 
fested in the long-range interaction case by a nonzero 
value of C at intermediate time scale To. We can use C 
as an order parameter to characterize different dynamical 
behaviors. By continuously tuning the effective interac- 
tion range (indicated by the fitting parametera) with the 
beatnote frequency p, we find that (see Fig. |4j) prether- 
malization only takes place when a is smaller than a crit- 
ical value (oec ~ 1-3 for N = 16). For larger system size, 
the prethermalization-thcrmalization transition still oc- 
curs, but etc becomes smaller and the transition becomes 
sharper. For the particular case a — 0, the system has 
uniform coupling and its dynamics can be solved exactly. 
The exact solution shows that the system stays in the 
prethermal state forever with C = % — 1. 

Mechanism of prethermalization and dynamical tran- 
sition . We now give a physical explanation to why 
prethermalization and different dynamical behaviors oc- 
cur in this model. The distinctive short time dynamics 
of (crf(t)) (Fig. [3^i, b) can be explained by examining 
the energy spectrum of Hxy in the single spin excita- 
tion subspace (shown in Fig. |5^i). In the short-range 
interaction case, the energy spectrum is close to linear. 
This is because Hxy can be roughly approximated with 
only neighboring interaction, which is then identical to a 



tight-binding fcrmionic model 

H tb = 2^(J M+i cjc m + Bc\Ci) 

i 

and shows "quantum mirror" behavior |291 l3U| , resulting 
in a near dispersion-free spin wave propagation until non- 
linearity sets in. On the other hand, the energy spectrum 
for the long range interaction case is highly non-linear, 
so the dynamics of spin excitation is strongly dispersive. 

The cause for existence of prethermal stage in the 
long time dynamics, however, is much more complicated. 
Naively, the spin flip-flop matrix J^- varies smoothly 
among sites for any a € (0,3), so the spin excitation 
should continuously diffuse from one end of the chain to 
the whole chain, and is not expected to get trapped some- 
where in the middle for a long time. The two stage dy- 
namics indicates that there are two different time scales 
weaved in the Hamiltonian, which is not at all obvious by 
looking at J^-. We note, however, that the time dynamics 
of any physical observable is simply given by 

(A(t)) =J2Pmn(0)A nm e i ^- E ^ h , 

where p mn (0) is the initial state's density matrix element 
in energy basis. So different time scales of dynamics can 
be unraveled through mapping of eigenenergy differences 
{E m — E n }, as done in Fig. [5] In the short-range inter- 
action case (Fig. ^p), all {E m — E n } are continuously 
distributed from J to 100 Jo, so a single-stage relaxation 
is expected after T ~ To = 10/ Jo- In the long-range in- 
teraction case (Fig. [5J3) , most {E m — E n } still fall into the 
range of 1 — 100 Jo , but there is a striking separate branch 
gaped at much lower rate (~ 10 -6 Jo). This branch actu- 
ally corresponds to near-degenerate pairs ({E2k~ •Eafc-i}) 
of eigenenergy (Fig. [5^) that make up the first half en- 
ergy spectrum, and the number of these pairs scale up 
with system size N. The appearance of near-degenerate 
pairs {E^k — i^2fc— 1}) in eigenspectrum of our model 
seems to be due to a combined effect of long-range in- 
teraction and inhomogeneous lattice spacing. If we put 
the ions into a ring (or flat-bottomed) trap so that the 
ions are equally spaced, we find that there is no sepa- 
rate branch in {E m — E n } plot (figure not shown) and 
hence no prethermalization behavior even with the same 
long-range interaction. 

Figure |5ji shows that the thermal values can be well 
predicted by the diagonal ensemble (DE) , defined as 

Pde — Pmn(0)5 mn , 

In large N limit, under certain conditions prescribed as 
eigenstate thermalization hypothesis [3H5], the diagonal 
ensemble prediction will match the canonical ensemble 
for thermalized state in classical statistical physics. 

As prethermalization is due to a different scale of 
eigenenergy difference in the Hamiltonian, to predict the 
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Figure 5: (a) Energy spectrum in the single spin excitation 
subspace for A = 16 spins with short and long range interac- 
tion, (b-c) Scatter plot of eigenenergy differences {E m — E n } 
for N = 256 spins with (b) short range interaction (a = 2.4) 
and (c) long range interaction (a = 0.74). (d) Comparison of 
prethermal values of af (blue circle, taken at t — To = 500ms) 
with PDE prediction (blue solid line) and thermal values (red 
square, taken at t — 5 x 10 3 s) with DE prediction (red dashed 
line) based on long range interaction pattern Jy shown in Fig. 
[2Jd. See in the text for the definition of PDE and DE. 



prethermal values here, we can define the partial diagonal 
ensemble (PDE): 

(0)6 I > l/Tb 

I Pmn (0) K-f„|«l/T 

where {u m } are the eigenenergies. We find that the PDE 
can well predict prethermal values of local observables 
of, as shown in Fig. [HJi. Roughly speaking, the prether- 
malization time scale is determined by the average level 
spacing (~ 10/ Jo), and the thermalization time scale is 
determined by the minimum level spacing (~ 10 4 / Jo as 
in Fig. |Ji). 

The dynamical transition can be associated with 
breaking of the lattice inversion (parity) symmetry, rem- 
iniscent of symmetry breaking in equilibrium phase tran- 
sitions. Our Hamiltonian Hxy is symmetric under the 
space inversion around z = (Fig. [lj, but we start from 
an initial state that does not have this symmetry. The 
thermal state, with no memory of initial state, is de- 
scribed by the diagonal ensemble Pde and restores this 
symmetry as C = 0. However, the prethermal state does 
not restore the Hamiltonian symmetry due to its non- 
zero C value, which indicates that some "memory" of 
initial state is preserved in prethermal state. The inter- 



mediate time scale To for observation of the prethermal 
state gives a microscopic interpretation why this state 
can break the parity symmetry: one cannot distinguish 
the near-degenerate pairs of eigenstates in the energy 
spectrum, so linear combinations within each pair are 
allowed. Since the two eigenstates of the pair have ei- 
ther even or odd parity, their linear combinations can 
break the parity symmetry. The dynamical phase dia- 
gram shown in Fig. [4} also has the hint of two non- 
analytic points: one is where C becomes non-zero, rep- 
resenting the appearance of prethermalization, and the 
other is where C approaches — 1 (a — > 0) , representing 
the disappearance of thermalization. 

Discussion of experimental detection: The transverse 
field Ising Hamiltonian Eq. [I] has already been exper- 
imentally simulated in Ref. [19| for N — 16 ions, with 
demonstrated highly-efficient in-situ measurement of spin 
polarization (of) and spin correlation (afcr|). The XY 
Hamiltonian (Eq. [2| can thus be readily obtained by 
tuning up the effective transverse magnetic field. The 
non-equilibrium initial state preparation requires focused 
laser beam, but is relatively easy due to large ion spac- 
ing near the ends. The laser power and trap frequencies 
used for generating interaction pattern Jy as shown in 
Fig. [2] are within current experimental reach |17l IT5] . 
The observation of prethermalization and the dynam- 
ical transition shown in Fig. [3] & [4] only requires the 
spin decoherence time longer than To = 10/ Jo = 500ms, 
and coherence time up to 2.5s has been experimentally 
achieved using the hyperfinc qubit of Yb + ions |31| . But 
the second stage of thermalization for long-range inter- 
action case will take much longer time and is beyond 
current experimental reach, similar to the experiment on 
prethermalization with cold atoms |14j . 

In summary, we have proposed a novel scheme to ob- 
serve the peculiar prethermalization phenomenon and 
dynamical transitions in the experimental system of 
trapped ion quantum simulator. The required conditions 
fit well with the current experimental technology. We 
provide an explanation of the mechanism of prethermal- 
ization and dynamical transition in our proposed model, 
which is connected with some unique feature of this ex- 
perimental system. 
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